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I ' We present the leading unpolarized and single spin asymmetries in inclusive two-hadron 

production in electron-positron annihilation at the Z pole. The azimuthal dependence in 

' the unpolarized differential cross section of almost back-to-back hadrons is a leading cos(2(/!>) 

^ ' asymmetry, which arises solely due to the intrinsic transverse momenta of the quarks. An 

OO ■ extensive discussion on how to measure this asymmetry and the accompanying time-reversal 

OO ' odd fragmentation functions is given. A simple estimate indicates that the asymmetry could 

nI ' be of the order of a percent. 

I \ 

^ ■ 13.65.-Hi, 13.85.Ni, 13.87.Fh, 13.88.-he 

; 

] Recently, we have presented the results of the complete tree-level calculation of inclusive two-hadron produc- 
• tion in electron-positron annihilation via one photon up to subleading order in 1/Q 0], where the scale Q is 
Oh! defined by the (timelike) photon momentum q (with = q^) and given by Q = y/s. The quantity Q had to 
^ ' be much larger than characteristic hadronic scales, but - being interested in effects at subleading order - we 
^ ' considered energies only well below the threshold for the production of Z bosons. 

f~| , In this article we extend those results to electron-positron annihilation into a Z boson, such that the results can 
' be used to analyze LEP-I data. We will neglect contributions from photon exchange and 'y-Z interference terms, 
which are known to be numerically irrelevant on the Z pole. Only leading order (1/Q)° effects are discussed, 
since for Q>Mz the power corrections of order 1/Q are expected to be completely negligible. Furthermore, 
^ , we will focus on tree level, i.e., order (as)*^- A rich structure nevertheless arises when taking into account the 
intrinsic transverse momentum of the quarks and, possibly, polarization of the detected hadrons in the final 
state. By accounting for intrinsic transverse momentum effects we extend the results in the analysis of Chen et 
al. where no azimuthal asymmetries arising from transverse momenta have been considered. 

For details of the calculation and the formalism we refer to We shortly repeat the essentials. We consider 
the process e~ -t- — s- hadrons, where the two leptons (with momentum / for the and I' for the e^, 
respectively) annihilate into a Z boson with momentum q = l + l', which is timelike with q^ = Q^. Denoting the 
momentum of the two outgoing hadrons by Ph [h = 1, 2) we use invariants Zh = 2Ph ■ q/Q'^ . We will consider 
the case where the two hadronic momenta Pi and P2 do not belong to the same jet (i.e.. Pi ■ P2 is of order 
Q^)- In principle, the momenta can also be considered as the jet momenta themselves, but then effects due to 
intrinsic transverse momentum will be absent. We will treat the production of hadrons of which the spin states 
are characterized by a spin vector Su {h = 1, 2), satisfying Ph ■ Sh — and — 1 < S'^ < 0. In this way we can 
treat the case of unpolarized final states or final state hadrons with spin-0 and spin- 1/2. In the present article 
we will disregard the polarization of hadron two (summation over spins). The final states which have to be 
identified and analyzed for the effects we discuss are simpler than the ones investigated by Artru and Collins 
[^, who proposed to measure azimuthal correlations in four-hadron production. 

The cross section (including a factor 1/2 from averaging over incoming polarizations) for two-particle inclusive 
e~^e~ annihilation is given by 
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PP pP c/cr(e+e~) q,2 

^1 ^2 ^ r ynf^-^ (I) 

d^P, (PP2 4 ((02 _ Af|)2 + r|Af|) Q2^M<''''' ' 

with Qfu, = e2/(167rsin^ 6w cos^ Svf) and the hchcity-conserving lepton tensor (neglecting the lepton masses and 
polarization) is given by 

L^.u[l. n = (5y' + 9a^) [^IX + ^U'^ - Q^9^^A - {^9v9\)'^^ e^.pJ'l'^, (2) 

where gy, 9^ denote the vector and axial- vector couplings of the Z boson to the leptons, respectively, and the 
hadron tensor is given by 



d'Px_ 
with 

H^APx;PiSi;P2S2) - {0\J,,mPx; PiSi; P2S2){Px; PiSi; P2S2\M0m, (4) 

where a summation over spins of the unobserved owi-state X is understood. 

In order to expand the lepton and hadron tensors in terms of independent Lorentz structures, it is convenient 
to work with vectors orthogonal to q. A normalized timelike vector t is defined by the boson momentum q and 
a normalized spacelike vector z is defined by P'' = — {P ■ q / q^) q'^'- for one of the outgoing momenta, say P2, 

0^ 

''^Q^ (5) 

(1 ~ P'^ 

Vectors orthogonal to t and z are obtained with help of the tensors 

5'^'' =g^'' + (7) 

e_L = -e T-pZa - 7^ re"^ P2pqa- (oj 

yP2 ■ q) 

Since we have chosen hadron two to define the longitudinal direction, the momentum of hadron one can be used 
to express the directions orthogonal to t and z. We define the normalized vector = P^^/\Pi±\, with P[Y = 
9^" Piu, and the second orthogonal direction is given by e'^h^, (see Fig. |l|). We use boldface vectors to denote 
the two-dimensional Euclidean part of a four- vector, such that Pi± ■ Pi± = —Pi± ■ Pij_. 




In the calculation of the hadron tensor it will be convenient to define lightlike directions using the hadronic 
(or jet) momenta. The momenta can then be parametrized (remember that Pi • P2 is of order Q^) using 
dimensionless lightlike vectors and 7i_ satisfying — r?_ — and rij^ ■ n_ = 1, 
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(9) 
(10) 

(11) 



with = —Q^. We have neglected hadron mass terms and considering the case of two back-to-back jets we 
have <C Q^. We will use the notation = p-riz^ for a generic momentum p. As momentum P2 defines the 
vector z^, 

Vectors transverse to n+ and n_ one obtains using the tensors 



(12) 



9t = 9 -nl'^nj, 



uv _ Ltupa 



(13) 
(14) 



where the brackets around the indices indicate symmetrization. The lightlike directions can easily be expressed 
int,z and a perpendicular vector, 
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V2 



[i^ + z"] , 

t> _ |M + 2 



(15) 
(16) 



showing that the differences between g'^ and gl^'^ are of order 1/Q. We will see however that taking transverse 
momentum into account does not automatically lead to suppression. 

To leading order the expression for the hadron tensor, including quarks and antiquarks, is 



W^"" = 3 y dp-dk+(fpj,(fkj, S^iprp + kj,- Qj,) Tr (A(p)F^A(fc)F'' 
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where = gvl^ + 9Al5l^ is the Z boson-quark vertex. We have omitted flavor indices and summation. The 
correlation functions A and A are given by 



S e**'-^ (0 1 (x) I Pi , ^1 ; X) (Pi , ^1 ; X I (0) 1 0) , 
'^x e-'P-^ (0|V^^-(0)|P2, ^2; ^)(P2, ^2; ^|^.(a;)|0) 



(18) 
(19) 



and the quark momentum k (and similarly for p) and the polarization vector (from now on we omit S'2) are 
decomposed as follows: 



z 



(20) 
(21) 



To leading order in l/Q one has that z — zi. The (partly integrated) correlation function A is parametrized as: 
i j dk+ A{PuSi;k) 
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(22) 
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(23) 



We parametrize the antiquark correlation function A in the same way, except that the distribution functions 
are overhned and the obvious replacements of momenta are done. 

The functions Di, . . . in Eq. ( p^ ) and Gil,Git, ■ ■ ■ in Gis,... are called fragmentation functions. One 
wants to express the fragmentation functions in terms of the hadron momentum, hence, the arguments of the 
fragmentation functions are chosen to be the lightconc (momentum) fraction z = /k~ of the produced hadron 
with respect to the fragmenting quark and kip = —zkrp^ which is the transverse momentum of the hadron in 
a frame where the quark has no transverse momentum. In order to switch from quark to hadron transverse 
momentum a Lorentz transformation leaving k~ and unchanged needs to be performed. The fragmentation 
functions are real and in fact, depend on z and k'j? only. 

We note that after integration over kj, several functions disappear. In the case of Tr(Aicr'~75) a specific 
combination remains, namely Hi = Hit + {k\/2Mi) H^. 

The choice of factors in the definition of fragmentation functions is such that J dz cPk'rp Di{z,k'rp) = N^, 
where is the number of produced hadrons. 

Note that the decay probability for an unpolarized quark with non-zero transverse momentum can lead to 
a transverse polarization in the production of spin-1/2 particles. This polarization is orthogonal to the quark 
transverse momentum and the probability is given by the function D^j,. In the same way, oppositely transversely 
polarized quarks with non-zero transverse momentum can produce unpolarized hadrons or spinless particles, 
with different probabilities. In other words: there can be a preference for one or the other transverse polarization 
direction of the quark (aligned or opposite relative to its transverse momentum) to fragment into an unpolarized 
hadron. This difference is described by the function H^. It is the one appearing in the so-called Collins effect 
1^, which predicts a single transverse spin asymmetry in for instance semi- inclusive DIS, and arises due to 
intrinsic transverse momentum. 

The functions D^rp and H^ are what are generally called 'time-reversal odd' functions. For a discussion on 
the meaning of this, we refer to Ref. |Q] and earlier references ; here we only remark that it does not signal 
a violation of time-reversal invariance of the theory, but rather the presence of final state interactions. 

The cross sections are obtained from the hadron tensor after contraction with the lepton tensor 



where {fiv} indicates symmetrization of indices, [^v] indicates antisymmetrization. The fraction y is defined to 
he y = P2 ■ I/P2 ■ q ~ l~ /q^ , which in the lepton center of mass frame equals y = (1 -I- cos02)/2, where 62 is the 
angle of hadron two with respect to the momentum of the incoming leptons. 

Azimuthal angles inside the perpendicular plane are defined with respect to l'^, defined to be the normalized 
perpendicular part of the lepton momentum l, l^ = l'^/{Q\/y{l ~ y))- 



{l-2y + 2y2) 5/^- + 4y(l - 2/)z^z 




(25Ua) (1 - 2y) 4" - 2« v/y(l - V) hpe^z 



(24) 



\a±_ \ cos (/), 



(25) 




\a±\ sin 0, 



(26) 



for a generic vector a. 



The vector and axial couplings to the Z boson are given by: 




(27) 
(28) 
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4. = (qv^ ~ 9a^) > j = £ or u,d,s 



4 = ^9Wa- (29) 
As well, we will use the following kinematical factors: 

^(y)= Q-y + y') = (i + cos2 02)/4, 

Biy)^y{l-y) sin^^sM 

C(y) = (1-2?;) -cos 02. (30) 

We obtain in leading order in 1/Q and the following expression for the cross section in case of unpolarized 
(or spinless) final state hadrons: 

dndzM^q^ ~^ iQ'- MlY + r|Af| '''' \ ^^^^ ^ ^3^(2/)) [d,u,\ 



-cos{2(f>i) 44 B{y) T < ^ ^ ^ '- - \ ^i^^i 



(31) 



M1M2 

where dfi = 2dyd4>^ and 0' gives the orientation of We use the convolution notation 

J'[D''W]= J d^kj. d^PT 5'^{pT + kj,-qj,)D''{zi, zl fc|) D° (z2 , ^aPr ) ■ (32) 

The angle 4>i is the azimuthal angle of h (see Fig. Q). In order to deconvolute these expressions we can define 
weighted cross sections 

= / ^ d\, w '"i^Trl^'^''^ (33) 

^ y 27r dQdzidz2d^qj, ' 

where W — W{Qt, 0i, 02, 0Si, '/'S2)- The subscript A denotes the polarization in the final state for hadron one, 
with as possibilities unpolarized (O, including the case of summation over spin), longitudinally polarized (L) or 
transversely polarized (T) . We postpone the discussion of the additional structures and information accessible 
by measuring the polarization of one of the final state hadrons to the end of this letter. 

Even without determining polarization of a final state hadron a subtle test of our understanding of spin transfer 
mechanisms in perturbative QCD can be done. The information on the production of a transversely polarized 
quark-antiquark pair, which subsequently fragment into unpolarized (or spinless) hadrons with probabilities 
depending on the orientation of the (anti)quark's spin vector relative to its transverse momentum, is contained 
in the cos(2(/)i) azimuthal asymmetryQ To access this information we utilize the weighted cross sections 

^^^o = _ ,^2 E (c^? A(y) - i 4 cgC(y)) I??(.i) I5?(z2), (34) 

where the fe^-moments for a generic fragmentation function F are defined by 



^This asymmetry is not to be confused with the cos{2(f>) asymmetry found by Berger JtI, which is suppressed. 
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Fiz,,zfkl). 



(36) 



We now like to focus on the weighted cross section defined in Eq. ( p5| ) and discuss its possible measurement. 
In order to be able to observe the cos(2(/)) dependence one must look at two jet events in unpolarized electron- 
positron scattering. In each jet one identifies a fast hadron with momentum fractions zi and Z2 respectively. 
One of the hadrons (say two) together with the leptons determines the lepton scattering plane as is indicated 
in Fig. |l|. In the lepton CM system hadron two determines the z-direction with respect to which the azimuthal 
angles are measured. One needs in particular the azimuthal angle (jfi of the other hadron (one) as well as its 
transverse momentum which determines Qt — \Pi±\/zi. The cos(2(/)) angular dependence then can be 

analyzed by calculating the weighted cross section of Eq. (|35| ) . 

For an order of magnitude estimate, we consider the situation of the produced hadrons being a tt"*" and a 7r~. 
Furthermore, we assume £'"""^^(2) = £'f"^^^(z) {Df^"^ (z) = Z?"^'^ (z), respectively) and neglect unfavored 
fragmentation functions like _D''^^^(z) etc.; and similar for the time-reversal odd functions. The equalities for 
the Di functions seem quite safe on grounds of isospin and charge conjugation, the same assumptions might be 
non-trivial for the functions. As a consequence of these assumptions the fragmentation functions can be 
taken outside the flavor summation, and we obtain 



Qt 



AM1M2 



Di(zi) £'1(2:2) 



where 



Fiy) 



(37) 



(38) 



This factor is shown in in Fig. ^ as a function of the center of mass angle 6*2 (we use sin^ 9w — 0.2315 [|j). 
At an angle close to 90° we observe the largest effect. In order to get an estimate of the true asymmetry 
at the level of count rates, one should compare Eq. ( |35| ) with the weighted cross section ((5|,/4MiM2)q. To 
estimate the ratio of those two quantities, we use as argued in reference ^ , for the ratio of the fragmentation 

functions H^'"^\zi) / Di{zi) = 0(1), although this is likely an optimistic estimate. From the average transverse 
momentum squared of produced pions in one jet, for which we take 0.5 (GeV/c)^ [|lO|, one obtains an estimate 
for the average transverse momentum of pions in jet one with respect to a given pion in jet two. This leads 
at zi = Z2 = 1/2 to (^Q"^ / AM^') ^ « 50 (1)q and consequently to an estimate at the percent level for the 
ratio (((5|,/4M^) cos(2(/)i))q / ^(3|,/4M^)^. Such an azimuthal dependence in the unpolarized cross section, 
however, may be detectable in present-day electron-positron scattering experiments. 



F(e,) 



30 60 90 120 150 180 "2 




FIG. 2. Factor defined in Eq. (|38|) depending on the center of mass angle 62. 



The situation where hadron two is taken to be a jet, which in this back-to-back jet situation is equivalent 
to analyzing the azimuthal structure of hadrons inside a jet, is obtained by considering D1{z2) = (5(1 — Z2) 
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and ^ (^2) = 0. This gives the familiar result for J dz2 and it gives zero for the Z2-integrated cos2(/) 
azimuthal asymmetry. 

The experimental determination of the polarization of (one of) the final state hadrons offers further oppor- 
tunities to reveal the hadronic structure in terms of spin-dependent fragmentation functions. We assume in 
the following that the spin vector of hadron one, i.e. Si, is known (reconstructed), having in mind the example 
of a produced A and its self- analyzing properties. We observe a rich structure of angular dependences due to 
polarization. 

Again, weighted cross sections are the appropriate means to separate out specific functions. For instance, the 
weighted cross section 

(S -('^^ +^^.))^ - (Q._g.^;r|,,| E c[c^. Biy) H^izi)Ht^nz.) (39) 

picks out the term which is the closest analogue to the Collins effect [|] in semi-inclusive lepton-hadron scattering 
prf . We note that a confirmation of the cos(20) asymmetry, also implies a confirmation of the Collins effect. 
A complete list of weighted cross sections at leading order is given in Table |. 

TABLE I. Weighted cross sections for Si / 0, 5*2 = 



W A (W)^ ■ [3ai q7 ((Q2 - M|)2 + r|M|)] 



1 


L 


iQ^/AMiMi) sm{2(t>i) 


L 


(Qt/Mi) sin((?!)i - (^s J 


T 


{Qr/Ah) sm{(t>i + (t)Si) 


T 


{Q%/&M'iM2) sin(30i - 4>s^) 


T 


(Qt/Mi) cos{(j)i - <j,s,) 


T 



\SM Ea,, H Aiy) - 1 44 Ciy)) Dt^^^'^{z,)D't{z2) 

-\SM Eara ^1^2 B{y) Hti''^\zi)H^'''^\z,) 
I^itIE... {c^i^Aiy) -icUlC{y)) G«"(zi) :DJ (^2) 



In conclusion, we have presented the leading asymmetries in inclusive two-hadron production in electron- 
positron annihilation at the Z pole. We have investigated unpolarized and single spin asymmetries. We 
included the effects of intrinsic transverse momentum and in this sense our results are an extension of those of 
Ref. [^. The azimuthal dependence in the unpolarized differential cross section is a cos(2(/)) asymmetry, which 
arises solely due to the intrinsic transverse momenta of the quarks. An extensive discussion on how to measure 
this asymmetry and the accompanying time-reversal odd fragmentation functions is given. A simple estimate 
indicates that the asymmetry could be at the percent level, hence it can perhaps be observed in present-day 
electron-positron scattering experiments. In confirming the existence of this asymmetry one also confirms the 
Collins effect, without the need of a polarization measurement. 

Note added: In a preliminary study a similar correlation in back-to-back jets was already experimentally 
investigated. We find that it involves moments of the functions and , different from the ones in our 
correlation. In this study no significant result was found using the 1991 to 1994 LEP data. In this analysis 
three momenta in the final state need to be determined, namely besides two hadron momenta also the jet axis, 
and hence there are two azimuthal angles, (j) 8'iid </>', yielding a cos(0 -I- (j)') asymmetry. 
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